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SU(4) Einstein- Yang-Mills theory possesses sequences of static spherically symmetric globally 
regular and black hole solutions. Considering solutions with a purely magnetic gauge field, based 
on the 4-dimensional embedding of su(2) in su(4), these solutions are labelled by the node numbers 
(ni,ri2, 113) of the three gauge field functions ui, U2 and 113. We classify the various types of solutions 
in sequences and determine their limiting solutions. The limiting solutions of the sequences of neutral 
solutions carry charge, and the limiting solutions of the sequences of charged solutions carry higher 
charge. For sequences of black hole solutions with node structure (n,j,n) and (n,n,n), several 
distinct branches of solutions exist up to critical values of the horizon radius. We determine the 
critical behaviour for these sequences of solutions. We also consider SU(4) Einstein- Yang-Mills- 
dilaton theory and show that these sequences of solutions are analogous in most respects to the 
corresponding SU(4) Einstein- Yang-Mills sequences of solutions. 
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I. INTRODUCTION 



SU(2) Einstein- Yang-Mills (EYM) theory possesses a sequence of asymptotically flat static spherically symmetric 
solutions, which are globally regular pL Based on a purely magnetic gauge field ansatz, these solutions carry no 
charge. They are characterized by the node number n of the single magnetic gauge field function u. With increasing 
node number n, the sequence of neutral solutions tends to a limiting solution with magnetic charge P — 1, which 
consists of two parts, an inner oscillating part and an outer Reissner-Nordstr0m (RN) part 

Besides the globally regular solutions, SU(2) EYM theory possesses static spherically symmetric black hole solutions 
with non-trivial non-abelian gauge field configurations outside their regular event horizon || . For any value of the 
event horizon radius x^, there exists a corresponding sequence of black hole solutions, characterized by the node 
number n of the gauge field function. For fixed horizon radius Xh and increasing n, the sequence of neutral black 
hole solutions also tends to a limiting charged solution with magnetic charge P = 1 0]. For horizon radius Xu > 1, 
this limiting solution corresponds to an embedded RN solution with the same horizon radius xn. For < xn < 1, 
the limiting solution again consists of an oscillating part and a RN part . Interestingly, the only charged black hole 
solutions of SU(2) EYM theory are embedded RN solutions 0. 

In SU(3) EYM theory, besides embedded RN solutions Q charged static spherically symmetric black hole solutions 
exist, which possess non-trivial gauge field configurations outside their regular event horizon [HQ|. These black hole 
solutions have magnetic charge of norm P = \/3 residing in the su(3) Cartan subalgebra J^]. Like the RN solutions, 
the magnetically charged non-abelian black hole solutions exist only for horizon radius xn > P, where ih = P leads 
to extremal black hole solutions 0. These charged non-abelian SU(3) black hole solutions are characterized by the 
node number n of a single magnetic gauge field function, similar to the neutral non-abelian SU(2) black hole solutions. 
With increasing n and fixed horizon radius xn, the sequence of non-abelian black hole solutions with charge P = \/3 
tends to a solution with higher charge, P = 2 M. For horizon radius sh > 2, this limiting solution again corresponds 
to an embedded RN solution with the same horizon radius xh, whereas for ^/i < xh < 2 the limiting solution again 
consists of an oscillating part and a RN part . 

SU(3) EYM theory also possesses neutral globally regular and black hole solutions. The genuine neutral static 
spherically symmetric SU(3) EYM solutions are obtained by embedding the 3-dimcnsional representation of su(2) in 
s«(3) ]lO| , pT| ,[7[ . The purely magnetic gauge field ansatz then involves two functions, u\ and U2, and the solutions 
are labelled by the corresponding node numbers (711,712). For fixed horizon radius xh, a discrete set of regular 
neutral SU(3) EYM solutions is obtained [Q. These neutral solutions form sequences with node structure (77,77) 
and (7,7 + 77), with i fixed. With increasing 77, these sequences of solutions again tend to limiting solutions, carrying 
magnetic charge of norm P = 2 and P = V3, respectively. [RH]. For the (7,7 + 77) sequences, the limiting solution 
of black hole solutions with horizon radius xh > \/3 corresponds to a charged non-abelian SU(3) black hole solution 
with i nodes and with the same horizon radius xh, whereas for black hole solutions with horizon radius < xh < \/3 
and for globally regular solutions the limiting solution again consists of two parts . The (77, 77) sequences either 
represent genuine SU(3) solutions with u\ ^ ui or scaled SU(2) solutions with iti = U2- Their limiting solutions carry 
magnetic charge of norm P = 2, and represent embedded RN solutions for ih > 2. Unlike the (7, 7 + 77) sequences of 
black hole solutions, the genuine SU(3) (77, n) black hole solutions exist only for sufficiently small values of the horizon 
radius xh- At critical values of the horizon radius, they merge into the scaled SU(2) solutions. Therefore, of the two 
types of (77, 77) sequences only the scaled SU(2) solutions persist for large xh 0- 

Here we investigate the sequences of SU(4) EYM solutions. The genuine static spherically symmetric SU(4) solutions 
are obtained by embedding the 4-dimensional representation of su(2) in s?j(4) [jl(J[l4|. The purely magnetic gauge 
field ansatz then contains three functions, u\, ui and 773. The classification of the discrete set of regular solutions 
in sequences again involves their node numbers (rii, 772, 773) [fl2"|| . When some of the gauge field functions possess the 
same number of nodes, several solutions can exist for sufficiently small values of the horizon radius xh, leading to a 
complex critical behaviour of the solutions as functions of xh~- When all three gauge field functions are non-trivial, 
neutral solutions are obtained, whereas when one or more gauge field functions are identically zero, charged solutions 
arise. We determine the limiting solutions of the sequences of globally regular and black hole solutions and specify 
their magnetic charge. The extension to SU(N) EYM theory is straightforward. 

By including a dilaton field, all these solutions can be generalized to the solutions of Einstein- Yang-Mills-dilaton 
(EYMD) theory The classification of the solutions of EYMD theory and EYM theory is identical. However, 

the charged EYMD black hole solutions exist for arbitrarily small values of the horizon radius xh, and the limiting 
solutions of some EYMD sequences of solutions are embedded Einstein-Maxwell-dilaton (EMD) solutions 

In this paper we present in section II the SU(4) EYM action and the static spherically symmetric ansatz for the 
metric and the purely magnetic gauge field, based on the 4-dimensional embedding of su(2) in sw(4), and we determine 
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the boundary conditions for asymptotically flat globally regular and black hole solutions. In section III we discuss the 
charged SU(4) EYM black hole solutions. We classify them and construct numerically several sequences of charged 
black hole solutions. In particular, we demonstrate the exponential convergence of some of their properties. In section 
IV we classify the globally regular SU(4) EYM solutions in sequences. We consider several sequences in detail and 
present numerical results for them. The sequences of neutral SU(4) EYM black hole solutions are discussed in section 
V, where we consider in particular the complex pattern of bifurcations, occurring at critical values of the horizon 
radius. In section VI we briefly generalize our results to SU(4) EYMD theory and present some numerical examples. 
We give our conclusions in section VII. 

II. SU(4) EINSTEIN- YANG-MILLS EQUATIONS OF MOTION 
A. SU(4) Einstein- Yang-Mills action 

We consider ST/ (4) EYM theory with action 

S = S G + S M = J L G ^gd i x + J L M V zr gd i x , (1) 

La = , (2) 



16ttG 



gauge field Lagrangian 



field strength tensor 



L M = --Tt(F^F^) , (3) 



gauge field 



Fpv = d^A v - d u A ll - ie[A^ A v ] , (4) 



A, = \\ a A% , (5) 



and gauge coupling constant e. 

Variation of the action eq. ([l]) with respect to the metric g^ v leads to the Einstein equations 



with stress-energy tensor 



G>„ — R^ v — -g^R — %ttGT^ v , (6) 



dL 1 
V = g^L M - 2^ = 2Tr(F A1Q F^ 3 «' 3 - -g^F afi F^) , (7) 



and variation with respect to the gauge field A^ leads to the gauge field equations. 

B. Static spherically symmetric ansatze 

To construct static spherically symmetric solutions we employ Schwarzschild-like coordinates and adopt the static 
spherically symmetric metric 

ds 2 = g^dx^dx" = -A 2 JVdt 2 + Af^dr 2 + r 2 (dB 2 + sin 2 6d<p 2 ) , (8) 

with the metric functions A{r) and 

tf(r) = l-^. (9) 
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The static spherically symmetric ansatze for the gauge field of SU(4) EYM theory are based on the su(2) 
subalgebras of su(4). Here we consider the gauge field ansatz corresponding to the 4-dimensional embedding of su(2) 
in su{4) gij 



/ 3 



cos 



wi9 



M 2e 



uji<d cos 6d(j) 
~cosdd(f> W30 



with 



\ 




\ w 3 6 -3cos(9d0 / 



= + sin 



(10) 



(11) 



i. e. 



A a = A r = 



(12) 



The ansatz contains three gauge field functions 0Jj(r), j — 1,2,3, and leads to the field strength tensor components 



F r Q = d r Ag 



(13) 



and 



with 



F r< f, — d r A^ 



F H = (l/2e)diag(/ 1 ,...,/ 4 )sin0 , 



fj = - ^j-i +Sj , Sj = 2j - 5 , j = 1, . . . , 4 (w = cj 4 = 0) 



Another parametrization of the gauge field ansatz is given in 14 



(14) 
(15) 
(16) 



C. Field equations 



With the above ansatze we derive the set of EYM equations. The metric (|8j) yields for the tt and rr components 
of the Einstein equations 



Gtt — 



2m'A 2 N 



= 8irGT tt , 



(17) 



and 



G tt . 2 A' 



A 2 M 2 r A 



= 8TrGT rr , 



(18) 



where the prime indicates the derivative with respect to r. The static spherically symmetric ansatz for the fields ( fioj ) 



yields for the tt component of the stress-energy tensor T tt = —A 2 AfL 



M: 



Tt, 



T A 2 jV(AfG + V) , 



(19) 



and for the rr component 



T — 

-L rr — 



(Afg - v) , 



with 



3=1 



4 r 2 ^ - 3 

J'=l 



(20) 



(21) 
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We introduce the dimensionless coordinate 



the dimensionless mass function 



and the scaled gauge field functions |Tc|| 



with fj = -yju] - 7i-i«i-i + Sj (see eq. ©). 

For the gauge field functions we obtain the equations 

1 



(22) 



era 

(23) 



7i=i(4-j), ./ 1 • (24) 



y/li 

The above Einstein equations then yield for the metric functions the equations 

fi = MQ + V , (25) 

where the prime now indicates the derivative with respect to x and 



(AA/V,-) + ^a^(/i+i - = , (28) 
where the metric function A can be eliminated by means of eq. ( p6| ) to yield 

x 2 Nu 3 + 2{n - xV)u 3 + i(/ 3+1 - / 3 -)uj = . (29) 
We note the symmetry of the equations with respect to the transformation 

ui ->■ u 3 , u 3 -nil ■ (30) 



D. Boundary conditions 



We now consider the boundary conditions for the SU(4) EYM solutions. Globally regular solutions must satisfy 
boundary conditions at the origin. These are PQj 



MO) = o 



(31) 



and 



11,(0) =±1, j = 1,2,3. (32) 
For black hole solutions with a regular horizon with radius boundary conditions are imposed at xr. These read 

M(x H ) = or h(xk) = "j- , (33) 

and 



1 



M u i + ^2(/i+i-/jK' 







(34) 
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For charged black hole solutions the coefficient of u'j in eq. ([34|) may vanish, i. e. Af = 0. Then extremal black hole 
solutions are obtained, which satisfy 

(/i+i - />iU = ■ (35) 

At infinity black hole solutions and globally regular solutions satisfy the same set of boundary conditions. Asymp- 
totic flatness implies that the metric functions A and fj, both approach a constant at infinity. We here adopt 

^((5o) = f , (36) 

thus fixing the time coordinate. The mass of the solutions is given by /x(oo). Magnetically neutral solutions are 
obtained, when all gauge field functions are non-trivial and satisfy the boundary conditions 

tt,-(oo)=±l, 3 -1,2,3, (37) 

yielding an asymptotically vanishing field strength tensor component Fg^ = 0. No globally regular magnetically 
charged solutions exist. Magnetically charged black hole solutions are obtained, when one or more gauge field functions 
are identically zero, m = 0. Their non-vanishing gauge field functions uj ^ itj then approach constants Cj ^ ±1 at 
infinity §, 

Uj (oo) = Cj , (38) 

determined in section III. A. 



III. MAGNETICALLY CHARGED SU(4) EYM BLACK HOLE SOLUTIONS 

When one or more gauge field functions are identically zero, Uj(x) = 0, magnetically charged SU(4) EYM black hole 
solutions are obtained. These solutions are important in classifying the neutral SU(4) EYM solutions in sequences 
and in constructing their limiting solutions. Therefore we here consider these charged SU(4) EYM black hole solutions 
first. A general discussion of the magnetically charged SU(N) EYM solutions was given previously Isj. 



A. Classification of the magnetically charged SU(4) EYM black hole solutions 

We now classify the magnetically charged SU(4) EYM black hole solutions, obtained within the ansatz ( |T~0| ) ||. 
When one gauge field function is identically zero, uik = 0, the ansatz reduces to 



with Hk = cos ed ^ hk and 



A^dx" = 



A 



(4-fc) 



dx^ 



(39) 



(4-fe)l 



(*■) 



-fcl 



(40) 



(JV-fe) 



Here A^ (N = fc,4 — k) denotes the non-abelian static spherically symmetric ansatz for the su(N) subalgebra of 
su(4) (based on the A^-dimensional embedding of su(2)) and Tik represents the ansatz for the element hk of the Cartan 
subalgebra of su(4). 

The gauge field functions of the su(N) parts of the solutions satisfy the boundary conditions (eq. (|37|)) 



Ui(oo) = ±1 , 



(41) 



corresponding to neutral su(N) solutions Identifying the non- vanishing functions u>j of the su(A) ansatz with the 
corresponding functions uii of the non-abelian su(N) ansatze, 



G 



7l u- = nfju 2 j 



(42) 



then yields the asymptotic boundary conditions for the functions Uj 



Uj (oo) = c 3 = ±\jnh ] ■ (43) 

The charge of the solutions is carried by the Cartan subalgebra part of the gauge field. A solution based on the 
element of the Cartan subalgebra carries magnetic charge of norm P, 

P 2 = ^Tr h\ . (44) 

Expanding the element hk in terms of the basis {A„2_i | n — 2,3,4}, the charge can also be directly read off the 
expansion coefficients, 

N 

h k = Y / d^P n 2_ 1 X n2 _ 1 , (45) 

n=2 



where P n *_ x = yJ^H. The expansion coefficients, corresponding to cases 2a-c, are shown in Table 1. 

By applying these considerations again to the subalgebras su(N) of eq. (^9|), we obtain cases la-c of Table 1 , where 
two gauge field functions are identically zero. In the special case (case of Table I) where all gauge field functions 
are identically zero, an embedded RN solution is obtained with charge of norm P = a/IO- 

RN black hole solutions exist only for horizon radius xh > P, and the extremal RN solution has xh = P- As first 
observed for SU(3) EYM theory @], the same is true for charged non-abelian black hole solutions. Non-abelian black 
hole solutions with charge of norm P exist only for horizon radii sh > P- For extremal black hole solutions Af = 0, 
so the coefficient of u- in eq. ([54]) vanishes. This yields the boundary conditions || 

Uj(x n ) = Cj , (46) 

corresponding to Ui(x-n) = ±1. 



B. Numerical solutions 



All possible cases of magnetically charged SU(4) EYM black hole solutions are classified in Table 1. Since several 
cases are equivalent, fa ~ lb ~ Ic and 2a ~ 2c, three non-equivalent non-abelian cases remain, which are discussed 
below. 

The magnetically charged non-abelian static spherically symmetric solutions of SU(4) EYM theory are labelled by 
the nodes of their non-vanishing gauge field functions, i. e. by the nodes of the subsets of gauge field functions Ui, 
belonging to the subalgebras su(2) or su(3). These solutions form sequences, completely analogously to the sequences 
of neutral solutions in SU(2) and SU(3) EYM theory, which are classified in sequences by their node structure, n for 
SU(2) and + n) and (n, n) for SU(3). With increasing n, the sequences of neutral solutions converge to charged 
solutions. Similarly with increasing n, the sequences of charged solutions converge to limiting solutions with higher 
charge. When the limiting solution has charge of norm P, then one of the gauge field functions Ui becomes identically 
zero for x > P in the limiting solution. In particular, the limiting solutions of case la are embedded RN solutions 
with charge of norm P = for x > P. In contrast, the limiting solutions of cases 2a and 2b are non-abelian 
solutions, which carry charge of norm P = 3 and whose su(2) gauge field function has j nodes. Here only for j = 
embedded RN solutions with charge of norm P = 3 are obtained for x > P. 

In Table 2 we show the mass /i(oo) of the first few solutions of several sequences together with the mass of 
their limiting solutions. For every case of Table 1, two sequences are shown: one sequence with extremal horizon 
corresponding to dh = V9, xn — and ih = V& for cases la, 2a and 2b, respectively, and one sequence with 
horizon radius ih = \/T0. For cases 2a and 2b, the second non- vanishing function of the sequences shown has no 
node. Therefore, in all three cases, the limiting solutions are embedded RN solutions for x > P, where P = \/T0 for 
case la, and P = 3 for cases 2a and 2b. 

We demonstrate the convergence of the functions for the extremal solutions of case la in Figs. 1. Fig. la shows the 
gauge field function u\ for n = 1 — 7 with odd n. Since the limiting solution has charge of norm P = vlO, we need 
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to distinguish two regions, an inner region 3 < x < y/lO, beginning at the horizon, and an outer region x > VlO. 
Satisfying the boundary condition ui(xh) = 1, the function u\ becomes increasingly steep with increasing n in the 
inner region, reaching a limiting solution whose first node resides just before x = \/l0, the norm of the charge of the 
limiting solution. Beyond x = \/T0 the function approaches zero in an exponentially increasing region. In Fig. lb we 
show the metric function Af '■ With increasing node number n the function Af approaches a second zero, located at 
x = \/l0, the charge of the limiting solution. Beyond x = vlO the function Af approaches the metric function of the 
extremal RN solution with charge P — ylO. In Fig. lc we show the charge function P(x), 

P 2 (x) = 2x (n(po) - fi(x)) . (47) 

With increasing n the charge function P{x) tends to the charge of the limiting solution in an exponentially increasing 
region. 

All this is strongly reminiscent of the behaviour of the neutral globally regular SU(2) EYM solutions and their 
limiting solution Hfl] ■ The limiting solution of the neutral globally regular solutions has charge of norm P = 1 . It 
also consists of two regions, an inner region < x < 1, beginning at the origin, and an outer region x > 1. Satisfying 
the boundary condition u(0) = 1, the function u becomes increasingly steep with increasing n in the inner region, 
reaching a limiting solution whose first node resides just before x = 1, the charge of the limiting solution. Beyond 
1 = 1, the function u approaches zero in an exponentially increasing region. Similarly, with increasing node number 
n, the function Af appoaches a zero located at x = 1, the charge of the limiting solution, whereas beyond x = 1 it 
approaches the metric function of the extremal RN solution with charge P = 1. 

From the above similarity between the neutral globally regular solutions and their limiting solution with the charged 
extremal solutions and their limiting solution, we conjecture, that the charged extremal black hole solutions here 
play the role of the neutral globally regular solutions in determining the convergence properties of the functions. 
Considering neutral black hole solutions, we observed previously, that the location of the innermost node of the 
corresponding globally regular solutions with respect to the black hole event horizon is strongly indicative about the 
degree of convergence of the black hole solutions to their corresponding limiting solutions . 

As for the globally regular solutions, the convergence of the location of the innermost node of the charged extremal 
black hole solutions is exponential. This is seen in Fig. 2, where we present Az* as a function of the node number n, 
with Az* being defined as the deviation of the location of the innermost node of the n-th extremal black hole solution 
from the location of the innermost node of the limiting solution, 

A^ = 14-4,1. (48) 

Considering the logarithm of Az*, we observe, that its slope is the same for the solutions of case la and case 2b. In 
contrast, the slope is different for the solutions of case 2a. Interestingly, for the globally regular neutral SU(2) EYM 
solutions the logarithm of Az* also has the same slope as for the solutions of cases la and 2b. We conclude, that it 
is only the (sub) algebra of the non-abelian solutions, which determines the slope. 

Similarly, the deviation of the mass of the n-th extremal black hole solution from the mass of the corresponding 
limiting solution, 

A^i n = ^oo(oo) - ^„(oo) , (49) 

decays with the same exponent, when the extremal solutions correspond to the same non-abelian (sub)algebra, as 
seen in Fig. 3, where in addition to cases la and 2b also the extremal charged SU(5) EYM solutions (with su(2) 
subalgebra) and the globally regular SU(2) EYM solutions are shown. 



IV. REGULAR SU(4) EINSTEIN- YANG-MILLS SOLUTIONS 

Here we consider the static spherically symmetric globally regular solutions of SU(4) EYM theory, based on the 
ansatz (|l0|). These discrete |l2j globally regular SU(4) EYM solutions are magnetically neutral. They are analogous 
to the globally regular SU(2) and SU(3) EYM solutions, obtained previously The SU(4) EYM solutions can 

be labelled by the node numbers (m, 712,713) of the gauge field functions u\, U2 and 113. The solutions can then be 
classified into sequences. We first discuss several sequences which have embedded abelian limiting solutions with 
magnetic charge of norm P for x > P. Then we discuss the general case and classify the general solutions into 
sequences, presenting some numerical examples. 
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A. Sequences with embedded abelian limiting solutions 

Here we discuss those sequences of globally regular solutions, which for x > P tend to embedded abelian limiting 
solutions with magnetic charge of norm P. In these sequences a single index n characterizes the number of nodes of 
one, two or all three gauge field functions, while the remaining gauge field functions have zero nodes. With increasing 
node number n the gauge field functions with zero nodes tend towards some finite constant value for x > P, whereas 
the gauge field functions with n nodes tend to limiting functions, which vanish for x > P. Consequently the limiting 
solutions represent embedded charged abelian black holes for x > P, whose charge may be read off Table 1. 

1. (n, 0, 0) sequence 

In Figs. 4a-c we present the first few globally regular solutions of the sequence with node structure (n, 0, 0), n = 1 — 7, 
with odd n. The first member of the (1, 0, 0) sequence is the solution with the lowest mass of all SU(4) EYM solutions 
(within the ansatz @), (i(oc) = 1.723995 (see Table 3). 

The limiting solution of this sequence is classified by case 2c of Table 1. For n — *■ oo, the limiting solutions 
corresponding to case 2c of Table 1 are EYM solutions based on the su(3) subalgebra of sit (4), which carry magnetic 
charge of norm P = y6- For x > P the gauge field function u\ of the solutions of this sequence tends to zero, while 
the functions ui and 113 approach the SU(3) vacuum solution. Therefore for x > P in this case the limiting solution 
corresponds to an embedded abelian solution, an extremal RN solution with the same charge (of norm P = y/6). 

In Fig. 4a the gauge functions u\, u-i and 1*3 arc presented. For x > P the limiting gauge field functions are u\ = 0, 
1*2 = \/l/2 and u 3 = */2/3, whereas for x < P they are non-trivial functions. With increasing node number n, the 
location of the innermost node of the function u\ tends exponentially to a finite limiting value and an accumulation 
point appears at x = P. 

In Fig. 4b we present the metric function Af. For small node number n (n < 3) Af possesses only one minimum, while 
for larger n several local minima of Af develop. For n — ► 00, the global minimum 7V| m in = occurs at x{Af\ram) = P- 
It represents the event horizon of the embedded extremal RN solution, which corresponds to the limiting solution for 
x > P. 

In Fig. 4c we present the charge function P(x). For finite node number n all solutions are neutral, and P(x) decays 
asymptotically to zero. However, with increasing n the charge function P{x) tends to the charge of the limiting 
solution P = V6 in an exponentially increasing region. 

Since the equations of motion are symmetric under the interchange of u\ and U3, the solutions with node structure 
(0, 0, n) form an equivalent degenerate sequence. 

2. Other sequences 

We now consider the other sequences of globally regular solutions which tend to embedded abelian solutions with 
charge of norm P for x > P. Turning to the sequence with node structure (0, n, 0), the limiting solution of this 
sequence is classified by case 2b of Table 1. The limiting solutions corresponding to case 2b of Table 1 carry magnetic 
charge of norm P = VE. For x > P they represent EYM solutions based on the su(2) © su(2) subalgebra of su(4). 
Since the gauge field functions u\ and U3 of the solutions of this sequence approach the SU(2) vacuum solution, in 
this case the limiting solution corresponds to an embedded abelian solution for x > P, an extremal RN solution with 
charge P = y/8. The masses /i„(oo) of this sequence are shown in Table 3. In addition to this sequence with u\ = U3, 
a second sequence with node structure (0,n, 0) might exist with u% / U3. However, we do not find such solutions 
numerically. 

Next we consider the cases where two gauge field functions have a finite number of nodes and the third function has 
zero nodes. The first such case is the sequence with node structure (n, n, 0). The limiting solution of this sequence is 
classified by case lc of Table 1. The limiting solutions corresponding to case lc of Table 1 carry magnetic charge of 
norm P = a/9. For x > P they represent EYM solutions based on the su{2) subalgebra of sit(4). Since the gauge field 
function u 3 of the solutions of the (n, n, 0) sequence approaches the SU(2) vacuum solution, in this case the limiting 
solution corresponds to an embedded extremal RN solution with charge P = V9 for x > P. The masses fi n (oo) of 
this sequence are shown in Table 3 p9| . An equivalent degenerate sequence is (0,n, n). 

Considering solutions with node structure (n, 0, n), we find two types of solutions leading to two distinct sequences. 
In the solutions of the first type the gauge functions u\ and U3 are identical, whereas in the second type u\ ^113. 
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The limiting solutions of both sequences are identical and classified by case lb of Table 1. Analogously to case lc, 
the limiting solutions corresponding to case lb of Table 1 carry magnetic charge of norm P — V9 and represent EYM 
solutions for x > P, based on the su(2) subalgebra of su(4). Only for the (n, 0, n) sequences the limiting solution 
corresponds to an embedded extremal RN solution with P = y/9 for x > P. The masses /x n (oo) of both of these 
sequences are shown in Table 3. 

Last we consider the solutions with node structure (77,77,77). Here three distinct types of solutions exist. In the 
solutions of the first type all three gauge field functions are identical, u% = U2 = Us. These simply represent scaled 
SU(2) solutions, with scaling factor \/l0. In the solutions of the second type only two gauge field functions are 
identical, u\ =113, and in the solutions of the third type all three gauge field functions differ from each other. For the 
second and third type of solutions, for a given n several distinct non-degenerate solutions may exist. These solutions 
are discussed in detail in section V.B. The limiting solutions of the (n,n,n) sequences are classified by case of 
Table 1. Thus for n — > 00, the sequences tend to an embedded extremal RN solution with magnetic charge of norm 
P = ylO for x > P. The masses fi n (oo) for several (n, n, n) sequences are shown in Table 4. 

B. General sequences 

All sequences of neutral globally regular solutions not included in the previous subsection tend to non-abelian 
limiting solutions for x > P, which carry magnetic charge of norm P. The norm of the charge of the limiting solutions 
is the same as in the corresponding abelian cases. In this subsection we first present a classification scheme for the 
general globally regular solutions of SU(4) EYM theory (based on the ansatz (|lC|)) and then demonstrate the main 
features of these solutions by constructing several sequences numerically. 

1. Classification of the general solutions 

The general globally regular solutions are labelled by the node numbers (711,712,713) of the corresponding gauge 
field functions Ui,U2,Us. The case, where all node numbers are identical, has already been considered in the previous 
subsection. Assuming next, that the node numbers of any two gauge field functions (uj) are identical, we obtain four 
types of sequences, labelled by the fixed index k and the running index n. For fixed k, k = 1, 2, 3,..., and running 
71, 71 = 1, 2, 3,..., these four types of sequences are (k, k,n + k), (k, n + k, k), (k, n + k, n + k) and (n + k, k, n + k). 
For k = the sequences of the previous subsection are obtained. In this scheme any solution with two identical node 
numbers is included exactly once, apart from equivalent degenerate solutions and the cases, where solutions with both 
u\ 7^ 773 and ui = U3 exist. The limiting solutions of the sequences [k, n + k, n + k) and (n + k, k,n + k) correspond 
to cases la and lb, respectively, and the limiting solutions of the sequences (k, n + k,k) and (k, k,n + k) correspond 
to cases 2b and 2a, respectively (Table 1). 

When all three node numbers differ from each other, we obtain a similar classification, which now requires three 
distinct indices, /, k and n. For fixed indices I and k with I = 0, 1, 2,..., and k = 1, 2, 3,..., and running index 77, there 
are three types of sequences, (/, I + k, I + k + 71), (I + k, I, I + k + n) and (I, I + k + n, I + k). Any solution with three 
different node numbers is thus included exactly once in this scheme, apart from equivalent degenerate solutions. The 
limiting solutions of the (I, I + k, I + k + n) and (I + k, I, I + k + n) sequences correspond to case 2a of Table 1 and the 
limiting solutions of the (I, I + k + n, I + k) sequence corresponds to case 2b. 

It is straightforward to extend this classification to the globally regular solutions of SU(N) EYM theory. 

2. (n + 1, 1, 0) sequence and (n + 1, 0, 1) sequence 

We now discuss the properties of these general sequences. Their qualitative features are analogous to those of 
the sequences discussed in the previous subsection. In order to demonstrate this, we here present two numerically 
constructed sequences: (n + 1, 1, 0) and (77 + 1, 0, 1). 

In Fig. 5a we present the gauge field functions of the globally regular solutions of the sequence with node structure 
(77 + 1, 1, 0) for 77 = 1 — 5. The limiting solution of this sequence is an EYM solution based on the su(3) subalgebra of 
s?j(4), with magnetic charge of norm P — y/6, corresponding to case 2c of Table 1. For x > P the limiting solution is 
the extremal SU(3) EYM black hole solution with node structure (1,0). In Fig. 5b the metric function Af is shown. 
With increasing node number n the metric function M tends to zero at x = a/6, the charge of the limiting solution. 
The metric function Af has the same qualitative behaviour for the (77, + 1, 1, 0) sequence as for the (n, 0, 0) sequence. 
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The charge function P{x) is shown in Fig. 5c. Again, with increasing n it approaches P = v6 in an exponentially 
increasing region. 

The mass /i(oo) of the solutions of the (n+1, 1, 0) sequence again converges exponentially to the mass of the limiting 
solution. This is seen in Fig. 6, where Afj, n (eq. (ff9|)) is presented. Also shown in Fig. 6 is the mass of the solutions of 
the (n+1, 0, 1) sequence. For a given n, the solutions of the {n+1, 0, 1) sequence have a lower mass than the solutions 
of the (n + 1, 1,0) sequence, but both sequences tend to the same limiting solution. For comparison, also the mass 
of the solutions of the (n, 0,0) sequence is shown. All three sequences are based on the su(3) subalgebra of su(4), 
and the function lnA/i n has the same slope for all three sequences. This is analogous to the case of the solutions 
based on the su(2) subalgebra of sit(4), discussed in section III.B. The masses of the solutions of the (n + 1, 1, 0) and 
(ri + 1, 0, 1) sequences and of their limiting solution are shown in Table 3. 



V. NEUTRAL SU(4) EINSTEIN- YANG-MILLS BLACK HOLES 



The SU(4) EYM black hole solutions are obtained ana logously to the globally regular solutions, but with boundary 
conditions imposed at the event horizon (see section II D ). As for the SU(3) EYM black hole solutions Jll[], in general 
sequences of neutral SU(4) EYM black hole solutions exist for all values of the horizon radius. However, for the SU(3) 
EYM solutions with node structure (n,n) there are two distinct types of solutions, the scaled SU(2) solutions with 
U\ = ui and the genuine SU(3) solutions with u\ 7^ u%. The genuine SU(3) solutions exist only up to a critical value 
of the horizon radius Xjf n , where they merge into the scaled SU(2) solutions. Analogously, in SU(4) EYM theory 
critical values of the horizon radius occur for the various types of solutions with node structure (n,j, n) and (n, n, n). 
Their complex critical behaviour is analyzed below fll9[|. 



A. General sequences 



To all of the globally regular solutions discussed above, the corresponding neutral black hole solutions exist at least 
for sufficiently small values of the horizon radius xr-- For fixed horizon radius xr, the sequences of neutral black hole 
solutions tend to limiting solutions with charge of norm P. For xr < P the sequences of black hole solutions tend 
to a non-trivial limiting solution for xr < x < P and to a charged non-abelian black hole solution or an embedded 
extremal RN solution for x > P. For xr > P the sequences of black hole solutions tend to a charged non-abelian 
black hole solution or an embedded RN solution with charge P and the same value of the horizon radius. Since the 
other qualitative features of the black hole solutions are completely analogous to those of the corresponding globally 
regular solutions, we do not discuss them here further. Instead we turn to the solutions with node structure (n, j, n) 
and (n,n, n), whose complex critical behaviour shows interesting novel features. 



B. Bifurcations 



We recall that for the globally regular solutions with node structure (n, 0, n), two types of solutions exist, one type 
with u\ — U3 and a second type with u\ 7^ 7Z3, whereas for the globally regular solutions with node structure (n, n, n), 
three types of solutions exist, one type with u\ — U2 — 113, representing scaled SU(2) solutions, a second type with 
u\ = 773 7^ U2, and a third type with u\ 7^ 7/3 7^ u-i- 

Similarly, for the black hole solutions with node structure (n, 0, n), there are two distinct types of black hole solutions 
for sufficiently small values of the event horizon. At a critical value of the horizon radius x§ n , the u\ 7^ U3 solutions 
with node number n merge into the u\ = 113 solutions. Beyond this critical value only the u\ = 773 solutions persist. 
With increasing node number n the critical value of the horizon radius x^ n decreases. For the black hole solutions 
with n = 1, 3 and 5 the critical values of the horizon radius are x^ x = 1.972, = 1.244 and x^ r 5 = 1.180 (see 
Table 5). 

In Fig. 7a the value of the gauge field functions at the horizon of the black hole solutions with node structure 
(n, 0,n) is shown for n = 1, 3 and 5, together with the corresponding three critical values of the horizon radius. For 
ri = 3 and 5 the ui = U3 solutions are only shown in the inset, in order not to cover part of the corresponding u\ 7^ M3 
solutions. For the u\ 7^ 113 solutions, we obtain two degenerate branches of solutions. They are degenerate because 
of the symmetry with respect to the interchange of u\ and 7*3, i. e. if in a given solution Tti(xn) lies on the upper 
branch and 773 (xh) lies on the lower branch, then in the degenerate solution 771 (xh) and 7x3 (xr-) are interchanged. In 
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Fig. 7b we present the corresponding masses. The mass of the u\ 7^ U3 solutions is always larger than the mass of the 
u\ = 113 solutions. In contrast, the temperature is lower for the u\ 7^ u% solutions, as seen in Fig. 7c. 

We now consider the black hole solutions with node structure (n, n,n), whose critical behaviour is more complex. 
As for the globally regular solutions, three distinct types of black hole solutions exist for sufficiently small values of the 
horizon radius. There is one solution of the first type with u\ = 1*3 = U2, we find two non-degenerate solutions of the 
second type with u\ = 113 7^ U2, and we find up to four non-degenerate solutions of the third type with u\ 7^ 113 7^ 1*2- 
At the largest critical value of the horizon radius x^' n , the two non-degenerate branches of ui = 1*3 7^ 112 solutions 
merge into each other. Notably, this critical value of the horizon radius i?j is slightly beyond the value of the horizon 
radius, where the branches intersect the branch of scaled SU(2) solutions. Beyond i^ 1 only the scaled SU(2) solutions 
persist. At the second largest critical value of the horizon radius xS' , two degenerate branches of ui ^ M3 ^ 112 
solutions merge into the branch of scaled SU(2) solutions. At the third critical value ar^' another two degenerate 
branches of u\ =/= U3 7^ 112 solutions merge into the lower branch of the u\ = 113 7^ U2 solutions. At the smallest fourth 
critical value i^ 4 two non-degenerate branches of Ui 7^ M3 7^ 1*2 solutions merge into each other, being neither close 
to the branch of scaled SU(2) solutions or one of the u\ — 7/3 7^ U2 branches. They cease to exist beyond . 

We illustrate this complex critical behaviour for the black hole solutions with node structure (1,1,1) in Figs. 8. The 
value of the gauge field functions at the horizon radius xn is shown in Fig. 8a for all solutions found, together with 
the four critical values of the horizon radius. The scaled SU(2) solutions, having u\ = U3 = U2, exist for all values of 
the horizon radius. We see two distinct non-degenerate branches of u\ — U3 7^ 112 solutions, up to the critical value of 
the horizon radius — 4.158. Of the six branches of «i 7^ 113 7^ U2 solutions two degenerate branches of solutions 
merge into the branch of scaled SU(2) solutions at the critical value = 2.714, another two degenerate branches 
merge into the lower branch of the u\ = U3 7^ 112 solutions at the critical value x^f = 1.025, as seen in the upper 
right inset of Fig. 8a, and two non-degenerate branches merge into each other at x^f = 0.445, as seen in the lower 
left inset of Fig. 8a. In Fig. 8b we present the mass of all branches of solutions as a function of the horizon radius. Of 
all seven distinct branches of solutions, the scaled SU(2) solutions have the lowest mass, and first one and then the 
other of the two branches of u\ = 1*3 7^ 112 solutions has the highest mass, except very close to x'S \ ■ The intersection 
of the two u% = 1*3 7^ U2 branches is seen clearly in Fig. 8c, where the temperature of all branches is shown. 

The critical behaviour of the black hole solutions with node structure (3, 3, 3) is shown in Fig. 9. The two non- 
degenerate branches of u\ — U3 7^ U2 solutions cease to exist at the critical value x^^ — 2.836, as seen in the upper 
right inset, two degenerate branches of u\ 7^ 113 7^ 112 solutions merge into the branch of scaled SU(2) solutions at the 
critical value s^ 2 = 2.090, as seen in the lower left inset, and two degenerate branches of u\ 7^ 1*3 7^ U2 solutions 
merge into the lower branch of u\ = 113 7^ 112 solutions at x c ^ — 0.474. The critical values of the horizon radius and 
the corresponding masses of the (n, n, n) sequences are presented in Table 4 for n = 1 — 5 pt|. 



VI. SU(4) EINSTEIN- YANG-MILLS-DILATON SOLUTIONS 

We here briefly consider the solutions of SU(4) EYMD theory. In EYMD theory the dilaton coupling constant 7 
represents an additional parameter. For 7 = EYM theory is recovered, whereas for 7 = 1 contact with the low- 
energy effective action of string theory is made. Most of the qualitative features of the static spherically symmetric 
EYMD solutions agree with those of the static spherically symmetric EYM solutions. This was observed previously 
for the solutions of SU(2) and SU(3) EYMD theory and it also holds for the solutions of SU(4) EYMD 

theory. In particular, the solutions of SU(4) EYMD theory can be classified in the same way as the SU(4) EYM 
solutions. Qualitative differences concern the charged black hole solutions and the limiting solutions of the sequences 
of solutions. As in EMD theory, where only the extremal solution with ih = has a naked singularity at the origin, 
in EYMD theory charged black hole solutions with a regular horizon exist for any value of the horizon radius xh > 0. 
The limiting solutions of the sequences of EYMD black hole solutions are EYMD and embedded EMD black hole 
solutions with the same horizon radius and the same dilaton coupling constant, which carry the same charge P as 
the limiting solutions of the corresponding sequences of EYM solutions, as demonstrated in detail for SU(3) EYMD 
theory Q. The limiting solutions of the sequences of globally regular EYMD solutions are extremal EYMD and 
embedded EMD solutions |l|fl|. 
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A. Action and Equations 



We briefly discuss the action, the equations of motion and the boundary conditions for the static spherically 
symmetric solutions of SU(4) EYMD theory. The SU(4) EYMD action is 

S = S G + S M = J L G ^d 4 x + J LZV=gd 4 x , (50) 

with Lq given in eq. (|J) and 

Lm = -\d,&&>* - exp^iTr^FH , (51) 

with field strength tensor F^ u given in eq. (Q), and dilaton coupling constant k. 

In EYMD theory we proceed analogously to EYM theory. We employ Schwarzschild-like coordinates and adopt the 
static spherically symmetric metric (|^). The ansatz for the static spherically symmetric gauge field is given by 
(|io|). The corresponding ansatz for the static spherically symmetric dilaton field is <f> = <f>(r). With these ansatze we 
obtain for the matter Lagrangian 

L n = _^_^a w+V)t (52) 



with Af, Q and V given in (|) and @. 

Changing to the dimensionless coordinate x given in ( p2[ ) and the dimensionless mass function fj,(x) given in (|2^), 
we further introduce the dimensionless dilaton field <b. 



4> = V47G$ , (53) 

and the dimensionless coupling constant 7, 



7 = n/V^G . (54) 

The choice 7 = 1 corresponds to string theory, while 4 + n dimensional Kaluza-Klein theory has 7 2 = (2 + n)/n . 
From the Einstein equations we obtain the field equations for n(x) and A(x), 

ft = ljAfx 2 <j>' 2 + exp(2 7 0) (AfG + V) , (55) 

"4 2/1 2 1/2 



A =~ ( 2 X ^ + ex P( 2 ^^)^J > ( 56 ) 

where the prime now indicates the derivative with respect to x. 

For the gauge field functions Uj(x), j = 1 — 3, and for the dilaton function cj>(x) we obtain the field equations 



exp(270) 
2x~ 2 

{ANx 2 4/)' = 2 1 A exp(2 1( f)) (Afg + V) . (58) 



(exp(2 7 0)^)' = ^flA(f j+1 - fj) Uj , (57) 



The metric and gauge field functions satisfy the same boundary conditions as in EYM theory, given in section 
II. D. These are supplemented with the boundary conditions for the dilaton function. At infinity the dilaton function 
satisfies 

0(oo) = , (59) 
for the globally regular solutions the boundary condition at the origin is 

0'(O) = , (60) 
and for the black hole solutions regularity at the horizon requires 

w= 2 7 exp(270) {Afg + v) (61) 
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B. Numerical solutions 



In SU(4) EYMD theory, sequences of neutral static spherically symmetric globally regular solutions can be con- 
structed for arbitrary dilaton coupling constant 7. Parameterized by the node numbers (rii, 112, TI3) of the gauge field 
functions (7x1,7x2,7x3), these sequences converge to charged static spherically symmetric extremal EYMD solutions, 
based on the subalgebras of su(4) as classified in Table 1 for EYM theory. 

We present some examples of globally regular solutions in the upper part of Table 5. Shown is the mass of the first 
few solutions (with odd n) of the sequences with node structure (n, 0, n) with ui = 1*3 and txi 7^ 1x3 for EYMD theory 
with 7=1 and for EYM theory (7 = 0). Also shown is the mass of the corresponding limiting solutions, representing 
an embedded extremal EMD solution for EYMD theory and an embedded extremal RN solution for x > P for EYM 
theory. 

Analogously, sequences of neutral static spherically symmetric black hole solutions can be constructed for arbitrary 
dilaton coupling constant 7. These sequences converge to the corresponding charged static spherically symmetric 
EYMD and embedded EMD black hole solutions. The temperature shown in the upper part of Table 5 for the 
globally regular solutions should be interpreted as the temperature of the corresponding black hole solutions in the 
limit ih - * 0. 

The lower part of Table 5 shows the mass and temperature of the black hole solutions of the (n, 0, n) sequences at 
the critical values of the horizon radius x^ a , where the u\ 7^ 7x3 solutions merge into the txi = 1x3 solutions, for EYMD 
theory with 7=1 and for EYM theory (7 = 0). Also shown are the critical values themselves and the value of the 
function u\ at x^ n . For 7 = 1 the critical values are smaller than for 7 = 0. This is in agreement with the critical 
values of the SU(3) EYMD black hole solutions 0, which with increasing 7 first decrease and then increase again. 

The charged EYMD black hole solutions also follow the classification scheme of the corresponding charged EYM 
black hole solutions. Since these were discussed at length for SU(3) EYMD theory in 0, we do not present further 
results here. 

VII. CONCLUSIONS 

We have investigated the various sequences of globally regular and black hole solutions of SU(4) EYM theory, 
based on the static spherically symmetric ansatz for the purely magnetic gauge field, obtained by embedding the 
4-dimensional representation of su(2) in su(4) ]l0| . We have classified the discrete set of solutions in sequences by 
means of the node numbers (rxi, 77.2, 113) of their gauge field functions u\, 1x2 and 7x3. We have determined the limiting 
solutions of the sequences of solutions, and we have constructed numerous sequences and their limiting solutions 
numerically. In particular, we have shown, that the exponential convergence of the mass of the globally regular and 
extremal black hole solutions depends on the non-abelian (sub)algebra of the solutions. 

The SU(4) EYM black hole solutions with node structure (n, j, n) and (n, n, n, ) show a complex critical behaviour. 
For such sequences several branches of black hole solutions exist up to critical values of the horizon radius. For 
sufficiently small values of the horizon radius two types of solutions with node structure (n, 0, n) exist. The u\ 7^ 7x3 
solutions merge into the u\ = 7x3 solutions at a critical value of the horizon radius x^ n . Beyond xg' n only the u\ = 113 
solutions persist. Again for sufficiently small values of the horizon radius there are three types of solutions with 
node structure (n,rx, n), scaled SU(2) solutions with u\ = 7x2 =1x3, solutions with u\ = 7x3 7^ 2x2, and solutions with 
u\ 7^ M3 7^ 7X2. The two non-degenerate branches of u\ = 7x3 7^ 1x2 solutions merge into each other and cease to 
exist beyond x^' n . At x^'^, two degenerate branches of txi 7^ 7x3 7^ 1x2 solutions merge into the branch of scaled 
SU(2) solutions. Another two degenerate branches of txi 7^ M3 7^ M2 solutions merge into the lower branch of the 
xti = 7x3 7^ U2 solutions at ij n 3 , and two non-degenerate branches of txi 7^ 7x3 7^ 7*2 solutions merge into each other at 

For SU(4) EYM theory, the classification of the solutions in sequences represents a generalization with respect to 
SU(3) EYM theory, since three node numbers are needed for the classification instead of two. It is straightforward 
to generalize the classification of the static spherically symmetric globally regular and black hole solutions further to 
SU(N) EYM theory, where (N — 1) node numbers are needed for the classification since (N — 1) gauge field functions 
are present. For the charged SU(N) EYM black hole solutions, the general classification is given in ||. An existence 
proof for regular SU(N) EYM solutions has been proposed recently jy|], which follows closely the existence proof for 
regular SU(2) EYM solutions 

For EYMD theory the same classification of the solutions holds. In EYM theory charged black hole solutions exist 
only for horizon radius ih > P, where P is the norm of the charge, whereas in EYMD theory charged black hole 
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solutions exist for any value of the horizon radius ih > as in EMD theory. The limiting solutions of the sequences 
of EYMD black hole solutions are charged EYMD and embedded EMD black hole solutions with the same horizon 
radius, whereas the limiting solutions of sequences of EYM black hole solutions are charged EYM and embedded RN 
black hole solutions with the same horizon radius for ih > P ■ The norm of the charge of the limiting solutions is 
determined by their su{A) subalgebras and is the same for EYM and EYMD theory. 

Recently new types of solutions were found in SU(2) EYM and EYMD theory. These are static globally regular 
and black hole solutions with only axial symmetry |2^,|23| and non-static non- rotational black hole solutions p4] , p5| . 
Such solutions as well as solutions with only discrete symmetries p(| are also expected in SU(N) EYM and EYMD 
theory with N > 2. 
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3 P 8 



Table 1 

The classification of the charged black hole solutions of SU(4) EYM theory is presented. Shown are the non- vanishing 
gauge field functions (denoted by Uj) and the identically vanishing gauge field functions (denoted by zero), the norm 
squared of the charge of the black hole solutions, P 2 , and the subalgebra of the solutions including the non-abelian 
subalgebra and the coefficients* (in the given basis) of the corresponding element of the Cartan subalgebra. 
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/i(oo) 






la 


2a 


2b 


u/xh 


^- 1 


1 ^- 1 


' ^- 


1 


3.12998 


2.83874 


2.96587 


2 


3.15685 


2.95799 


2.99427 


3 


3.16139 


2.98951 


2.99906 


4 


3.16213 


2.99741 


2.99985 


5 


3.16225 


2.99936 


2.99997 


7 


3.16228 


2.99996 


3.00000 


oo 


3.16228 


3.0 


3.0 




I H = Vio 


1 


3.13343 


2.89417 


2.98031 


2 


3.15864 


2.98734 


2.98734 


3 


3.16191 


3.00255 


3.00412 


4 


3.16224 


3.00406 


3.00416 


5 


3.16226 


3.00415 


3.00416 


7 


3.16228 


3.00416 


3.00416 


oo 


3.16228 


3.00416 


3.00416 



Table 2 

The dimensionless mass /x(oo) of the first few SU(4) EYM black hole solutions of the sequences corresponding to 
case la of Table 1, case 2a with node structure (n,0) and case 2b with node structure n and is presented for their 
respective extremal horizons and for the horizon radius xh = VTO. For each sequence the corresponding limiting 
value of the mass is shown in the last row (denoted by oo). 
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Moo) 




n 


0, n, 


nnO 


n,0,n 


n+1 1 


n+1, 0,1 




Ul / tl3 / 112 


1 


1.72400 


2.01903 


2.44237 


2.32822 




2.29105 


2.70316 


2.66330 


2 


2.21466 


2.57219 


2.86605 


2.82600 




2.82138 


2.79238 


2.78119 


3 


2.37577 


2.74891 


2.96878 


2.95668 




2.95588 


2.82453 


2.82109 


4 


2.42674 


2.80401 


2.99439 


2.98933 




2.98916 


2.83425 


2.83323 


5 


2.44252 


2.82096 


2.99879 


2.99738 




2.99734 


2.83740 


2.83707 


oo 


2.44949 


2.82843 


3.0 


3.0 




3.0 


2.83874 


2.83874 



Table 3 

The dimensionless mass n(oo) of the first few solutions of several globally regular SU(4) EYM sequences is presented. 
For each sequence the corresponding limiting value of the mass is shown in the last row (denoted by oo). 
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n 


1 


2 


3 


4 


5 


1 


cr , 1 
"^H .n 


4.1575 


2.9564 


2.8362 


2.8165 


2.8133 


/ircg(oo) (upper) 


2.7636 


3.0925 


3.1507 


3.1604 


3.1620 


/i rGg (00) (lower) 


2.7394 


3.0887 


3.1501 


3.1603 


3.1620 




3 1 874 


3.1322 






3.1621 


2 


cr 
X H .n 


2.7143 


2.1738 


2.0904 


2.0772 


2.0741 


/i rcg (oo) 


2.6734 


3.0780 


3.1484 


3.1600 


3.1619 


U cr,2(oo) 


2.8880 


3.1059 


3.1528 


3.1607 


3.1620 


3 


J H .n 


1.0246 


0.5595 


0.4741 






/ircg(oo) 


2.7436 


3.0889 


3.1502 






/i cr,3(oo) 


2.7816 


3.0925 


3.1507 






^rcg(oo) (Ul = U3 = W2) 


2.6204 


3.0717 


3.1475 


3.1599 


3.1619 



Table 4 

The bifurcation behaviour of the SU(4) EYM (n, n, n) sequences is presented for n = 1 — 5. The mass of the regular 
solutions is shown together with the corresponding values of the critical horizon and the mass of the black hole 
solutions at the critical horizon p0| . For the solutions with u\ = U3 =^ 1/2 (first set), there are two non-degenerate 
branches. For comparison, the mass of the scaled SU(2) solutions with u± = 113 = u-i is also presented. 
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(n,0,n) 


Xh — 




//(oo) 


T/Ts 


n 


7 = 





1 





1 


1 


Ml / U 3 


2.3282 


1.6269 


0.1031 


0.3445 


Ml = U 3 


2.2910 


1.5919 


0.1559 


0.4296 


3 


'"1 / ^3 


2.9567 


2.0886 


0.0075 


0.0954 




2.9559 


2.0875 


0.0094 


0.1104 


5 


Ml / U3 


2.9974 


2.1193 


0.0005 


0.0237 


Ml = M3 


2.9973 


2.1193 


0.0006 


0.0272 


oo 




3.0 


2.1215 


0.0 


0.0 



(n,0,n) 




//(oo) 


T/Ts 


ui (a&) 


n 7 = 





1 





1 





1 





1 


1 


1.9723 


1.7742 


2.5065 


2.0424 


0.3038 


0.5860 


0.8632 


0.8552 


3 


1.2440 


0.2927 


2.9633 


2.1058 


0.0154 


0.1411 


0.9063 


0.8812 


5 


1.1803 


0.0690 


2.9978 


2.1203 


0.0009 


0.0345 


0.9130 


0.8860 



Table 5 

Upper part: The dimcnsionless mass //(oo) of the first few u\ ^ u 3 and u\ — u 3 solutions with odd n is shown for 
the two globally regular sequences with node structure (n, 0, n) for SU(4) EYM theory and SU(4) EYMD theory with 
7 = 1. Also shown is the temperature T/Ts (Ts = (47t;z;h) _1 ), to be interpreted as the xu — * limit of the black hole 
solutions. The corresponding limiting values are shown in the last row (denoted by oo). 

Lower part: The critical values of the horizon radius x^ nl where the U\ ^ u 3 solutions with node structure (n,0,n) 
merge into the corresponding U\ — solutions are shown for SU(4) EYM theory and SU(4) EYMD theory with 
7 = 1. Also shown are the dimensionless mass //(oo), the temperature T/Ts an d ui(xg). 
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Fig. la: The gauge field function U\(x) is shown as a function of the dimensionless coordinate x for the charged 
SU(4) EYM black hole solutions of case la of Table 1 with extremal event horizon xh = 3 and node numbers n = 1, 
3, 5, 7. The inset illustrates the convergence of the innermost node to a value just below x = \/T0 (thin vertical line) 
in the limit n — > oo. The thin horizontal line represents the limiting function for x > VTO- 
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Fig. lb: Same as Fig. la for the metric function Af(x). The inset illustrates the emergence of a second zero at 
x = VlO in the limit n — > oo. 
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Fig. lc: Same as Fig. la for the charge function P(x). The thin horizontal line represents the norm of the charge 
of the limiting solution, P = a/10- 
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Fig. 2: The logarithm of the absolute deviation from the limiting solution Az* = — z\\ is shown as a function 
of the node number n for the location of the innermost node of the extremal SU(4) EYM black hole solutions of cases 
la, 2a and 2b of Table 1. 
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Fig. 3: The logarithm of the absolute deviation from the limiting solution Afi n = /^(oo) — /i„(oo) is shown as a 
function of the node number n for the mass of the extremal SU(4) EYM black hole solutions of cases la and 2b of 
Table 1, together with case la of ref. 0, representing extremal SU(5) EYM black hole solutions, as well as for the 
globally regular SU(2) EYM solutions. 
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Fig. 4a: The gauge field functions Uj{x), j = 1 — 3, are shown as functions of the dimcnsionless coordinate x for 
the globally regular SU(4) EYM solutions with node structure (n,0, 0) and node numbers n = 1, 3, 5, 7. The inset 
illustrates the convergence of the innermost node to a value just below x = \/& (thin vertical line) in the limit n — > oo. 
The thin horizontal line represents the limiting function for x > \/6. 
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Fig. 4b: Same as Fig. 4a for the metric function Af(x). The inset illustrates the emergence of a second zero at 
x = \/6 in the limit n — > oo. 
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Fig. 4c: Same as Fig. 4a for the charge function P(x). The thin horizontal line represents the norm of the charge 
of the limiting solution, P = ^/6. 
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Fig. 5a: The gauge field functions Uj(x), j = 1 — 3, are shown as functions of the dimcnsionless coordinate x for 
the globally regular SU(4) EYM solutions with node structure (n + 1, 1, 0) and node numbers n = 1, 3, 5. The thin 
horizontal line represents the limiting function for x > y/6. 



30 



etric Function (n+1,1,0 








6 



8 



log 10 (x) 



Fig. 5b: Same as Fig. 5a for the metric function Af(x) and node numbers n = 1 — 5. The inset illustrates the 
emergence of a second zero at x = \/6 (thin vertical line) in the limit n — > oo. 
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Fig. 5c: Same as Fig. 5a for the charge function P(x) and node numbers n = 1 — 5. The thin horizontal line 
represents the norm of the charge of the limiting solution, P = y/6. 
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Fig. 6: The logarithm of the absolute deviation from the limiting solution A/i„ = /^(oo) — /i„(oo) is shown as 
a function of the node number n for the mass of the sequences of globally regular SU(4) EYM solutions with node 
stucture (n, 0, 0), (n + 1, 1,0) and (n + 1, 0, 1). 
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Fig. 7a: The value of the gauge field functions U\ and u% at the horizon is shown as a function of the horizon radius 
Xh for the SU(4) EYM black hole solutions with node structure (n,0, n) and node numbers n = 1, 3, 5. For n = 1 
both m — u 3 and ui ^=113 solutions are shown. For n = 3 and n = 5 the ui = 113 solutions are only shown in the 
inset (Axh = 0.1, Am = 0.02). The critical values of the horizon radius xg n are also indicated. 
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Fig. 7b: Same as Fig. 7a for the dimensionless mass /u(oo). The inset (Axh = 0.7, A^i — 0.0003) only shows n = 5. 
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Fig. 7c: Same as Fig. 7a for the temperature T/T s (T s = 47r.T H ~ 1 ). The inset (Ax H = 0.5, AT/T S = 0.0003) only 
shows n = 5. 



3G 




Fig. 8a: The value of the gauge field functions U\, u 2 and U3 at the horizon is shown as a function of the horizon 
radius xn for the SU(4) EYM black hole solutions with node structure (1, 1, 1). The upper right inset (Axh = 0.03, 
Au = 0.004) shows the merging of the degenerate branches of u\ 7^ 113 7^ u 2 solutions (dotted lines) into the lower 
branch of u\ = u 3 7^ u 2 solutions at x^f. The lower left inset (Axh = 2 x 10~ 4 , Au = 9 x 10~ 6 ) shows two 
non-degenerate branches of u\ 7^ u 3 7^ u 3 solutions close to their critical value igj . 
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Fig. 8b: Same as Fig. 8a for the dimensionless mass n(oo). The inset (Axh = 0.45, A/i = 0.012) shows two 
non-degenerate branches of u\ ^ ^ ui solutions close to their critical value . 
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Temperature (1,1,1) 
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Fig. 8c: Same as Fig. 8a for the temperature T/Tg- The inset (Axjj = 0.3, AT/Tg = 0.02) shows two non-degenerate 
branches of «i ^ u 3 ^ u 2 solutions close to their critical value x^f. 
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Fig. 9: Same as Fig. 8a for the solution with node structure (3, 3, 3). The upper right inset (Axh = 0.05, Am = 0.07) 
shows the two non-degenerate branches of u\ = u a ^ u 2 solutions close to their critical value x^' 3 as well as the branch 
of scaled SU(2) solutions. The lower left inset (Axh = 0.1, Au = 0.05) shows two degenerate branches of ui ^ u 3 ^ u 2 
solutions close to their ciritical value well as the branch of scaled SU(2) solutions. 
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